Two approximate formulae for the binding energies in A— hypernuclei and 
light nuclei. 



Th.E.LioliosQ 

Department of Theoretical Physics, University of Thessaloniki,Thessaloniki 54006, Greece 
Abstract 

Two approximate formulae are given for the binding energies in A— hypernuclei and 
light nuclei by means of the (reduced) Poeschl- Teller and the Gaussian central potentials. 
CK ' Those easily programmable formulae combine the eigenvalues of the transformed Jacobi 
^ ■ eigenequation and an application of the hypervirial theorems. 

U' PACS: 21.80. +a,21.10.DT 
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A wealth of experimental results has been produced in an attempt to determine the 
binding energies of a A— hyperon in A— hypernuclei defined as: 



^, E={M,ore + MA)c'-MHc' (1) 

0^ ■ where Ma, Mcore and are the masses of the hypernucleus, its core and of the 

Q ! A— hyperon respectively. Those experiments range from the nuclear emulsions to the 

I strangeness exchange reactions [ffl 0, B| and up to the more recent associated production 



reaction. [|, |, g g 



O 

0> ' The importance of those energies stems from the fact that the A— particle, unlike other 

elementary particles (like pions or K- mesons), is a probe of the inner nuclear density and 
not of the nuclear surface. Moreover the binding effects of the A— hyperon influence the 
O I average properties of the nucleus such as the nuclear rms radius. An accurate theoretical 
^ I determination of the A— binding energies is, therefore, of paramount importance. In fact 
various methods have been employed using phenomenological A— nucleus potentials or self- 
consistent calculations. Although the derived theoretical formulae are satisfactory, more 
convenient methods are needed for a handy interpretation of the experimental results. That 
^ I need becomes more intense when it comes to the study of the mass number dependence of 
the energy quantities in A— hypernuclei and light nuclei. 
The reduced Poeschl- Teller (RPT) central potential : 

V{r) = ^ (2) 

cosh^ (^) 

and the Gaussian one: 



.,2 



V{r) = -Voe-T^ (3) 

can be used as approximations to the self consistent potential in the study of the behavior 
of a nucleon in light nuclei ([^)and have also been used extensively in Hypernuclear Physics 
(0] ~ EH)- fact the above potentials belong to the general class : 
'''emaihthcoholQphy sics.auth.gr 
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V{r) = -Vofi-) 



< r < oo 



(4) 



where Vq > is the potential depth, R > the potential radius and the "potential form 
factor" /, (/(O) = 1) , is an even analytic function of x = Namely, 



k=0 



2k 



where dk are the numbers: 



dk 



1 d' 



2k 



{2k)\dx^'' 



x=0 



k = 0,1,2,... 



di<0 



However, the corresponding Schroedinger eigenvalue problems: 
Poeschl- Teller: 



d^ 

dx^ 



. +g ^cosh ^{x) + Eni 



Unl = 



Unl (0) = Unl (oo) = 



Gaussian: 



(5) 



(6) 



(7) 



jf_ _ /(/ + !) 
dx'^ 



+ S-V^> Enl 



Unl (0) = Unl (oo) = 

where, x = j^, s = ( 2jvoR^ ) ^ Eni= 'S~^^,are not exactly solvable ( but for the 
s— states of the Poeschl- Teller one) 

By means of the Hypervirial theorems(HVT), a method used extensively in Theoretical 
Physics([|n| — [^), there were obtained approximate analytic expressions for the eigenval- 
ues of the potentials in question in the form of s— series. For instance for the RPT potential 
the eigenvalues are : 



r^HVT 

Enl = 



-1 + 2a„;S 



12 



12ali - Al{l + 1) + 3l + — [15 - Al{l + 1)] s^- 



945 



/(/ + 1) 12a^, - 4Z(/ + 1) + 3 



907200 



60 

{33280a^i/(/ + 1) 



■12816 [/(/ + 1)]' + 25 [440/(/ + 1) + 567]} s' - y^^^(^ + 1) {32720^,+ 



+Sali [1445 - 828/(/ + 1)] + [3 - + 1)] [1084/(/ + 1) - 2445]} / + .. 
where = 2n -|- / -|- 3 /2 

Moreover, there has been pointed out that the Jacobi eigenequation (JEQ) : 



(9) 
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sinh (a;) cosh ix 



+ 



Unl 



(10) 



Unl (0) = Unl (OO) = 



with 



'nl 



(Inl 
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can be used as an approximation to the Schroedinger eigenvalue problem for the Gaussian 
potential[^, thus obtaining very accurate approximate analytic eigenfunctions for the s— 
states of that potential [|l5[] . 

The similarity between ([TOD and (|^) is obvious . In fact if it wasn't for the centrifugal 
term of (^ they would have been identical ||13|. Note that, for the s-states {I = 0), series 
(|^) is summed up to the well known exact formula for the RPT potent ial[p!0[]: 



nO- 



'nO 



12) 



It is therefore plausible to expand (|TTD in a way similar to @ , that is with respect to s, 
and study the two expansions term by term. In fact, it turns out that if one disregards higher 
terms ( whose contribution is negligible anyway, especially for relatively heavy nuclei) the 
first terms of that HVT-series can be given by a very simple, easily programmable formula 
of a satisfactory accuracy. 

Therefore, after some algebra, the proposed approximate formula for the eigenvlaues of 
the RPT potential obtained is: 



RPT 
nl 



J L 

Snl + ^ 



1 16 

1 ttnis I a, 

5 105 



L] 



104 
945 



(13) 



where L = / (/ + 1) . 

In the same spirit one can use the JEQ in order to approximate the first terms of the 
corresponding HVT series for the eigenvalues of the Gaussian central potential. In that 
case the approximation is {very) satisfactory only for the {ground state) s-states . Namely 



e: 



Gauss 
nO 



J ^ 2 "^"^ 3 ^42 



(14) 



The kinetic energy can be readily derived by applying the Hellman-Feynmann theorem p7[ 
which in our case takes the form: 



< T > 



nl 



s dEni 
2 ds 



(15) 



To assess the accuracy of those formulae two tables have been set up which depict the 
energy eigenvalues for various A— hypernuclei calculated through ([T^) and ([T^) respectively, 
against those calculated 
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a) by a numerical integration of the corresponding Schroedinger equation. 

b) by means of a perturbation method |jl2[ 

c) experimentally ||28| 



The parameters of the potentials used here were obtained by fitting procedures to the 
available experimental data. 

In the tables we use the following notation: 

A =The mass number of the core nucleus of the hypernuclei 

n =The principal quantum number 

/ = The angular momentum quantum number 

E^^™ =The energy eigenvalues obtained by numerical integration 

-^n""' =The energy eigenvalues obtained by the present formulae 

E^^ =The energy eigenvalues( A— binding energies) obtained experimentally. 

E^f =The energy eigenvalues obtained through a perturbation method 

Hence, for the RPT potential we obtain the following table : 

TABLE I 

A n I E^^^^iRPT) 

28 
A 
31 
A 
51 
A 
51 
A 
8c 
A 
89 
A 



According to Table. I, the proposed formula for the RPT potential gives satisfactory 
results especially for the lower excited states of relatively heavy nuclei. Note that the per- 
turbation method gives no result for the 2p, 2d states (or higher ones with nodes). 

As for the Gaussian potential the corresponding table is: 

TABLE II 





A 


n 


/ 






T-iexp 


E'T (RPT) 


Si 


27 





1 


-6.86 


-6.63 


-7.0 ±1 


-7.02 


'S 


31 





1 


-7.76 


-7.56 


-8 ±0.5 


-7.92 


V 


50 





1 


-10.77 


-10.66 


-12± 1 


-11.05 


V 


50 





2 


-3.94 


-3.88 


-4.0 ± 1 


-4.11 




88 





1 


-14.07 


-14.01 


-15.5 ± 1.0 


-14.23 


'Y 


88 





2 


-7.69 


-7.64 


-9.5 ±1 


-7.90 




131 


1 


1 


-5.77 


-5.77 








131 


1 


2 


-1.67 


-1.74 









A 


n 


/ 


E^P' (Gauss) 


(Gauss) 


7-1 exp 


E^f (Gauss 




15 








-12.97 


-12.66 


-12.5 ±0.35 


12.39 




27 








-16.18 


-16.00 


-16.0 ±0.28 


-16.19 


32 c 
A "J 


31 








-16.87 


-16.71 


-17.5 ±0.5 


-16.88 


fv 


50 








-19.04 


-18.96 


-19.9 ± 1 


-19.05 


89v 
A ^ 


88 








-21.26 


-21.22 


-22.1 ± 1.6 


-21.27 




137 


1 





-11.02 


-11.36 




-11.05 



Obviously, for the Gaussian potential, the approximation is extremely good for the 
ground state while for the higher states there is room for improvement. The validity 
of the proposed formulae can, of course, be further assessed by their application to the 
experimental results [1-8], just as it was done in the above tables. 

A clear advantage of the present method, apart from its simplicity and its application 
to states not covered by the perturbation method in ref.[|l^ (see Table I), is its s and L de- 
pendence which provides a straightforward means of accuracy assessment. Namely Eq.(|l^) 
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shows that the lower the excited state the better the approximation , as the parameter L 
becomes smaller . 

On the other hand the s-dependence of the method isolates and clarifies the effects of 
the width and depth of the potential as well as of the mass number of the nucleus as shown 
in some recent studies[^, Thus, the deeper and wider the potential the smaller the 



parameter s and consequently the better the accuracy of the method. As for the mass 
number of the nucleus , the heavier the nucleus the smaller the parameter s and inevitably 
the better the accuracy. Moreover, in an effort to derive approximate wavefunctions of 
the potentials in question, one can equate the energy eigenvalues instead of the rms radii 
15], in order to use the parameter s as a fitting one. That approach is simpler but has 



not been tried yet. 



That s-dependence was also shown p9[ to enhance greatly the study of the mass number 



dependence of the energy quantities through the relation 

R = ro (A) ^1/3 (15) 

where R is the nuclear radius of a nucleus with mass number A, assuming a rigid core 
model. Note that another significant advantage is that the kinetic and the potential energy 
can be derived readily from and (|1^) by a simple differentiation with respect to s , 
through Eq.(|l5D. 



Of course there are remarkable applications of that s— dependent method in other fields 
such as in Atomic Physics]^ . 



We should underline that the present formulae could be also used both in the study of 
nucleon-nucleus interactions and whenever a handy fair approximation of the eigenvalues 
of the potentials considered is needed. 

Finally it would be very interesting, indeed, to apply the present method to the general 
class of potentials (§) and study the region of its validity which is expected to be very 
satisfactory for some particular form factors. 
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